Coupling of thermal and mass diffusion in regular binary thermal lattice-gases 
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We have constructed a regular binary thermal lattice-gas in which the thermal diffusion and mass 
diffusion are coupled and form two nonpropagating diffusive modes. The power spectrum is shown to 
be similar in structure as for the one in real fluids, in which the central peak becomes a combination 
of coupled entropy and concentration contributions. Our theoretical findings for the power spectra 
are confirmed by computer simulations performed on this model. 
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The power spectrum of light scattered by a binary so- 
lution is more complicated than that of a single compo- 
nent fluid. The central peak contains combined effects of 
entropy and concentration fluctuations ||] . The cross 
effects are well known in nonequilibrium thermodynam- 
ics as the Dufour effect and the Soret effect, and are 
caused by the coupling between heat-flow and diffusion. 
For example, heat may be transported by conduction but 
also by diffusion of the two components. In the simpli- 
fied two-dimensional model presented here this coupling 
phenomena can be analyzed in detail. 

Lattice gas automata (LGA) are plagued by many de- 
fects. Although some of these defects can be solved, the 
models are in general not suitable for modeling realistic 
fluids because they do not exhibit a fully realistic ther- 
modynamical behavior. They are, however, useful tools 
for understanding more fundamental problems in ther- 
modynamics related to discretization and testing con- 
cepts in statistical mechanics ||. In LGA the positions 
and velocities of point-like particles are discretized onto 
a lattice jfj. The dynamics consist of a cyclic process 
of propagating particles according to their velocity to a 
neighboring node followed by local collisions that typi- 
cally conserve mass and momentum. 

In order to recover the macroscopically isotropic 
Navier-Stokes equations the lattice in two-dimensions is 
usually chosen to be hexagonal. Apart from problems re- 
lated to this, early LGA were also plagued by unwanted 
spurious invariants || ||, [7]] . A satisfactory extension to 
include thermal properties in LGA was made by Grosfils, 
Boon, and Lallemand (GBL), who introduced a multiple 
speed model, defined on a two-dimensional hexagonal lat- 
tice [|). The model uses a velocity set consisting of a 
single rest particle and three rings, each containing six 



directions, with velocities of magnitude 1, y/3, and 2. 

Generalizations of LGA to mixtures have been used to 
study interfaces and phase transitions (^] . Some of these 
models use a passive label to distinguish the different 
species J10| . In other models the particles live on different 
lattices |ll], |l2j . However, these models are athermal and 
do not show the coupling phenomena we are interested in. 
A two species thermal model by using a passive label, is 
also not suitable, because mass and heat transport would 
completely decouple |fl3|| . 

Here we introduce, for the first time, a thermal bi- 
nary lattice-gas model capable of capturing the essence 
of a real mixture with respect to coupling of entropy and 
concentration fluctuations. This model allows us to cal- 
culate the coupling quantitatively and separate the differ- 
ent contributions, which in the continuous case in general 
only is possible in the low density limit. The two species 
of particles live on separate two-dimensional hexagonal 
lattices and are labeled red and blue. Independently they 
would behave as normal GBL models, but we allow the 
particles to interact during the collision phase, i.e. mo- 
mentum and/or energy can be transferred from one to 
the other lattice. The number of red and blue particles, 
however, is constant in every collision. 

The state of a node can be specified by a set of Boolean 
occupation numbers n.^, denoting the presence or ab- 
sence of a particle of type fi = {r, b} in velocity channel 
Ci, where i is a label running over all 19 velocities. Due 
to the Boolean nature of LGA, the ensemble average of 
the occupation numbers in equilibrium, is described 
by a Fermi-Dirac distribution 
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where a r , ctb, (3, and 7 are Lagrange multipliers. Here, 
(3 is the inverse temperature, a r and ah fulfill a chemical 
potential role, and 7 is a parameter conjugate to the flow 
velocity. For simplicity we will work in the overall zero 
momentum case by setting 7 = 0. 
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A collision outcome is chosen with equal probability 
amongst all members of an equivalence class, i.e. a group 
of states having the same red mass, blue mass, total mo- 
mentum, and total energy. In a regular binary mixture at 
most two particles, each of different type, can be in the 
same velocity state [ [l4[ . This differs from the color mix- 
ture where the particles of the original single specie model 
are given a color to distinguish between them and was 
analyzed in the case of the GBL color mixture (CGBL) 
jl3| . In that case the collision operator could be split into 
two separated steps: a GBL collision and an independent 
redistribution of the colors, which made the simulations 
a relative simple extension of normal GBL simulations, 
even though the number of different states in the CGBL 
model was 3 19 . Here we have a "true" 38- bits model 
leading to 2 38 states on which the collision operator has 
to act. Clearly a naive lookup-table strategy in order 
to simulate this system is only feasible if the table is 
constructed partially and stored temporarily during the 
simulations, due to the excessive memory requirements 
for storing the complete table. 

Regular binary mixtures, however, do allow for a con- 
venient solution, which is less efficient than storing the 
complete collision table, but has still a relative good per- 
formance (about a factor 2-5 slower than CGBL). Rather 
than storing a collision table based on the states, we make 
one based on the different classes. This allows us to gen- 
erate a set of outgoing classes for the two species and via 
a GBL-like process an outgoing state. If this is combined 
with the 48 symmetry operations (six rotations, two re- 
flections, red/blue symmetry, and particle/hole symme- 
try), we get a working algorithm that can be used on a 
computer with 256 MB of memory |l^] . 

The theoretical framework for thermal lattice-gases is 
well established |l6| , and is here extended to binary mix- 
tures. It is possible to solve the many-body dynamics 
of LGA by using the Boltzmann molecular chaos as- 
sumption. Although one would expect deviations for 
higher densities, it turns out that these deviations are 
very small, in fact for many purposes even within a few 
percent accuracy. Therefore, provided the fluctuations in 
the average occupation numbers are small, a Taylor ex- 
pansion of the collision term in the neighborhood of the 
equilibrium distribution is justified, yielding a linearized 
collision operator f2. Then, in first approximation the be- 
havior of the system can be obtained by analyzing these 
deviations in terms of eigenmodes. 

In the analysis of the behavior of fluctuations in LGA 
mixtures it is convenient to introduce the colored scalar 
inproduct (l| 

(A\B) = J2A(c itt )B( Cili )K ifi , (2) 
if/ 

where = /^(l — /i^) is the variance in the average oc- 
cupation number. For reasons of symmetry k is included 
in definition of the right vectors, i.e. \B)i^ — Hj U B(cj U ). 

Following the method of Resibois and Leener |18| we 
need to find the fc-dependent eigenfunctions and eigen- 
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FIG. 1: The diffusive transport values as a function of the 
density at reduced temperature 9 = 0.1 and fraction of red 
particles P r — 0.9. 

values of the single-time step Boltzmann propagator 

e- tk - c (l + n)\m) =e z(fc) hHfc)>- (3) 

Within this formulation the hydrodynamic modes are 
characterized by the fact that the eigenvalues z M should 
go to zero in the limit of small wavevectors. Therefore 
we can derive these fe-dependent eigenvectors and eigen- 
values by making a Taylor expansion in fc. The resulting 
eigenvalues to second order in fc are given by 

z±(fc) = ±ic s fc-rfc 2 , (4) 

z ± (k) = -vk 2 , (5) 
z s± (k) = -4fc 2 . (6) 

It is with these modes that the binary-mixture responds 
to deviations from thermal equilibrium. The first two 
eigenmodes describe sound propagating (in opposite di- 
rections parallel to fc) with T the sound damping coeffi- 
cient and c s the adiabatic sound speed. The third eigen- 
value describes a shear mode with v the shear viscosity 
coefficient. The last two eigenvalues represent purely dif- 
fusive, non-propagating processes, but are in general not 
directly related to a familiar transport coefficient. Rather 
they always appear in combination with each other and 
can be expressed in other transport properties by 

s° ± = \{x + -D)±y±Q 2 + (x-V)\ (7) 

where \ 1S the generalized thermal diffusivity, T> a prop- 
erty related to the mass diffusion, and Q can be con- 
sidered as a measure of the correlation between the two. 
Another interesting relation can be obtained by introduc- 
ing the ratio of the specific heats 7 

r = |(f + (7-i)x). (8) 
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In Fig. [y the diffusive transport properties are shown 
at fixed reduced temperature 8 = exp(— i/3). The modes 
Sj_ converge for low densities to the thermal diffusivity 
and mass diffusion, due to the decoupling of the fluctua- 
tions in entropy and energy. This behavior resembles the 
situation for binary solutions where a formula similar to 
(0) exists ^] , with the same decoupling in the very dilute 
limit. In our model this means that the ratio Q/(x ~ *D) 
vanishes. The value of Q, however, will in general remain 
small but finite due to the divergencies of the transport 
coefficients in the low density limit of LGA. 

The decoupling is also observed in the limit of a single 
specie, i.e. the fraction of blue particles goes to zero. 
In this limit the model reduces to a normal GBL model, 
albeit that the diffusion related property T> will remain 
finite. For reasons of symmetry s± also converge to \ 
and T> in the limit of equal red and blue density. In fact 
this special limit can be analyzed completely in a similar 
way as done for CGBL jl3|]. Finally decoupling appears 
in two other cases due to the duality of the model under 
interchanging particles and holes, these are the high den- 
sity limit and the limit in which one of the sublattices is 
almost completely filled. 

The decoupling arises from an effective equipartition 
in the model making the ratio of the average occupa- 
tion numbers fir/fa the same constant for each velocity 
channel i. As can be seen in Fig. [I] there seem to be also 
intermediate values at which decoupling appears. Al- 
though indeed the converge to x an d 2? this is not 
due to decoupling via equipartition but due to the can- 
cellation of terms. This property could obviously be used 
in the analysis, but the location of these points depends 
in a non-trivial way on the chosen system parameters. 

The Boltzmann approximation enables us to calculate 
the modes up to fairly large wavevectors, even within the 
generalized hydrodynamic regime. In the hydrodynam- 
ical regime of small wavevectors fc, small frequencies u>, 
the hydrodynamical modes are well separated from the 
kinetic modes that, due to their exponential decay, can be 
neglected. In combination with a Taylor-expansion, this 
allows one to derive a Landau-Placzek approximation of 
the spectral density S(k,uj) 

s(k,u) v 7-1 (, + x-t> \ 4fc 2 

S(k) ^ 7 V s° + - s°_ J to 2 + (s° ± k 2 ) 2 

x - i rfc 2 
+ ^7(w± Cs fc) 2 + (rfc 2 ) 2 

1 fc ^ c s k±uj 

+ 7 ~s ^ (LU±C s k) 2 + (Tk 2 ) 2 ' 

(9) 

where S(k) is the static structure factor. 

The spectrum contains an unshiftcd central peak and 
is formed by two Lorentzians due to the non-propagating 
processes characterized by s±. The shear mode does 
not contribute to the spectrum and the two propagating 
modes lead to the presence of the two frequency-shifted 
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FIG. 2: The central part of the spectrum in the Boltzmann 
approximation, the Landau-Placzek approximation, the s± 
contributions are shown separately. The system parameters 
are 9 = 0.1, p = 6.5, P r = 0.05, k x = 4 x 2tt/512. The 
wavevector k and frequency ui are given in reciprocal lattice 
and time respectively. 

Brillouin lines. The last two terms in Eq. (@) have a 
contribution orders of magnitude smaller than the am- 
plitude of the Lorentzians and induce a slight pulling of 
the peaks towards the central peak Q]. The symmetry 
of the different contributions is such that the ratio of 
the integrated contributions of the central peak and the 
Brillouin components is constant and given by 7 — 1. 



Fig. g illustrates the composition of the central peak. 
The simulation results overlap perfectly with the Boltz- 
mann and Landau-Placzek approximations. The contri- 
butions of the diffusive modes to the central peak are 
indicated separately and for comparison a least square 
fit based on Eq. (||) to the whole spectrum, including the 
Brillouin peaks outside the interval shown, is made based 
on approximating the central peak by a single Lorentz. 

As can be seen from the Landau-Placzek expression (||) 
one of these contributions will vanish in the limits where 
s± converge to \ and V and the central peak reduces to 

S™(k,Lj) -y-1 2 X fc 2 

S(k) 7 c^ 2 + (xfc 2 ) 2 ' 1 J 

even though the other transport coefficient remains fi- 
nite. As mentioned before, this occurs at a limited set 
of locations, such as the low/high density limit, single 
specie limits, and the limit of equal red and blue density. 
In addition to these general limits it is also found in the 
low/high temperature regions close to densities where all 
rings of particles with the same absolute velocity are com- 
pletely filled or empty |l3j . This reduction of the central 
peak to a single Lorentz in the low density limit is also 
found in the very dilute limit of binary solutions. But 
contrary to what is found here, it is usually the thermal 
diffusivity that disappears Q 0]. 
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FIG. 3: A diffusive spectrum based on the different pos- 
sible choices for the density difference. The curve labeled 
by * is corrected by subtracting the propagating part. The 
system parameters are = 0.05, p = 10.0, p r /p = 0.15, 
k x — 20 x 27r/512. The wavevector k and frequency ui are 
given in reciprocal lattice and time respectively. 



If s± are sufficiently different, one can separate the 
two contributions in the spectrum outside the limits men- 
tioned above (Fig. |J). This does not automatically im- 
ply that one could obtain the transport values of interest 
from the experimental spectrum only, because they can 
still differ significantly from x & n d T>. 

In the low density limit we can identify T> with the 
mass-diffusion of the GBL model. In the limit of a single 
specie this is no longer true if one goes to higher densi- 
ties, even though all other properties converge to their 
GBL values. The origin of this problem is found in the 
possibility to have more than one particle with the same 
velocity in this model. Continuous theory suggests that 
V is not the correct generalization of the diffusion [^J. 
This is consistent with the way in which this quantity 
arises in the present model Jig] . 



The proper generalization of the diffusion is not com- 
pletely trivial. There is some ambiguity in the choice one 
needs to make. In a continuous fluid one would consider 
the decay of a signal of the type 



1 r r b 



(11) 



where \R)i ii = S rfJ _ and \B)i ll = 5b^- A spectrum based 
on this normalized density difference, however, leads to 



Brillouin peaks as is shown in Fig. |^. This was already in- 
dicated earlier in an athermal binary mixture fl4|| , where 
it was proposed to use 



|diff) = 



\R) 

(R\R) 



\B) 



(BIB) 



(12) 



But in the thermal case this signal suffers from the same 
problem. In both cases one could in principle extract the 
propagating part, but it is more natural to reformulate 
the density of interest to 



|diff> 



\B) 



_M 

(p\R) (p\B) ' 



(13) 



with p = c 2 /2 being the energy per particle. This prob- 
lem originates from the lack of equipartition, but in gen- 
eral the differences between the various choices will re- 
main small []l5|| . 

In conclusion, we presented a thermal binary mixture, 
which exhibits spontaneous thermal and concentration 
fluctuations. The regular mixture model is much more 
involved than the color mixture, but it is closer to the dy- 
namics of real binary mixtures, as it is able to capture the 
coupling phenomena observed in binary solutions. The 
coupling between energy transport and diffusion results 
in two diffusive non-propagating modes and in a more 
complicated structure of the power spectrum. In gen- 
eral, however, these two "true" transport coefficients do 
not seem to correspond with a conventional transport 
coefficient, but always appear in combination with each 
other. 

The LGA spectrum is similar in structure as the one for 
the continuous case. In general it is impossible to sepa- 
rate the entropy and concentration contributions, and the 
central peak can not be described by a single Lorentzian. 
However in several limits, namely the low/high-density 
limit, the single specie limit, and the equal red/blue limit, 
the true modes s5j_ converge to generalized thermal diffu- 
sivity x an( i a mass-diffusion like T>. In other cases the 
different transport coefficients can be determined by us- 
ing the theoretical framework, that provides an accurate 
description as is demonstrated by the comparison of the 
Landau-Placzek expression with simulation results. 
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